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Abstract 

Large-scale structural patterns commonly occur in network models of complex sys- 
tems including a skewed node degree distribution and small-world topology. These 
patterns suggest common organizational constraints and similar functional conse- 
quences. Here, we investigate a structural pattern termed pathway proliferation. 
Previous research enumerating pathways that link species determined that as path- 
way length increases, the number of pathways tends to increase without bound. We 
hypothesize that this pathway proliferation influences the flow of energy, matter, 
and information in ecosystems. In this paper, we clarify the pathway proliferation 
concept, introduce a measure of the node-node proliferation rate, describe factors 
influencing the rate, and characterize it in 17 large empirical food-webs. During 
this investigation, we uncovered a modular organization within these systems. Over 
half of the food-webs were composed of one or more subgroups that were strongly 
connected internally, but weakly connected to the rest of the system. Further, these 
modules had distinct proliferation rates. We conclude that pathway proliferation in 
ecological networks reveals subgroups of species that will be functionally integrated 
through cyclic indirect effects. 
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1 Introduction and Motivation 



Large-scale structural patterns have been uncovered in network models of 
complex systems, suggesting the possibility of common organizational con- 
straints and similar functional consequences. Network models are mathemat- 
ical graphs composed of nodes and undirected edges or directed arcs that 
connect the nodes. For example, in a social network nodes might represent 
individuals in a community and the edges or links could r epresent a social re- 
lationship between the individuals such as collaboration ( Newman . l2001alfbh . 
In a r nodel of the World Wide Web, web pages are no des connected by hyper- 
links (|Albert et all Il999l : iBarabasi and Albertl . Il999f ) . Ecologists use network 
models in many ways, including to represent trophic rel ations in foo d -webs 



and more generally energy-matter flux i n ecosystems ( Cohen et al. . 199d 



Higashi and Burni . 1991 : Margalel, 19631 : Pimm . 1982| ). In these networks, 
species or functional groups form the node set while the presence of energy 
and matter transfers and transformations are represented by links. 



Trad i tional l y, complex systems have been modeled u sing random graphs (|Erdos and Renvi 
1959L Il96nl : FGardner and Ashbvl . Il97nl : iMavl . Il972[ ). However, ecologists have 
demonstrated that random graphs are inadequate models of ecological sys- 
tems; food-web and ecosy stem models o ften contain stru c tures not commoni 



found i n random graphs (ICohen et al. ', '1990'; 'DeAneelif^, '1975"; ' Lawleil . Il9^ 



Pimmi I1979L Il98j IPimm et al.rF991: Ulanowicz and Wolff. 19911) In food 



webs , these hypo thesized structures inclu de short food-chain lengths (jPimm and Lawton 
1977 : Post, 2002[ ) and little or no cycling ( Cohen et al. . 1990| ). In addition, sev- 
eral forms of modularity — hierarchic compartmentaliza tion into subsystems- 
have been hypot hesized for food-webs and ecosystems ("Allen an d Staril Il982l : 
Mav j_ I972I: Pimm. .l979i: Pirnm and Lawtort 1980.: Yodzis, 1982) . Furthermore , 
Ulanowicz and Wolfa (|l99lh demonstrated that random networks (based on 
the Poisson, uniform, Gaussian, negative exponential, and log-normal prob- 
ability distribution functions) failed to capture the distribution of connec- 
tions in real ecosystems. Ecologists hypothesize that these structural differ- 
ences exist because ecological systems are shaped and constrained by ther- 



modynamic laws and natural historv (lj0rgensen et al. . 19921 : Lawler . 19781 : 



Miiller and Leupeltl . Il998l : IWilliams and Martinez! . I2OOOI ) . 



Likewise, investigations of other types of complex systems have identified a 
number of disti nctive patterns common in complex systems not found in purely 
random graphs ((Albert and Barabasi , 20021 : Moreno and Jenningsl . 19381 : iNewman , 
2OO2L I2OO3I : [PricT Il965l ). For example, the distribution of node degree (i.e., 
the number of edges (links) incident to a node) is often skewed in models of 
complex systems, following an exponential distribution or a power-law distri- 
bution rather than the Poisson distribution of random graphs. The power- 
law distribution was found in the World Wide Web (^Barabasi and Albert . 
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19991 ), metabolic networks iTeong et all 120001) and some but not all food 
webs ( Dunne et al. . 2002bl : Montova and Sola 20021 ) . The power-law deg: 



ree 



distribution implies that there are a large number of nodes with ver y few con- 
necti ons, while a few nodes have a large number of connections (jBarabasi . 
2002[ ) . This topology tends to increase network robustness to ra ndom node or 



edge deletion, while making it more sensitive to targeted attacks ([Albert et al 



2000[ iDunne et al.'. '2002b'j^ The small-world pattern i s another commonly 



found topology (jWattSi . A99^ Watts and Strogat"3 . 19981 ). In small- world net- 
works, the degree of node clustering is larger and the maximum distance 
(where distance is the shortest path between two nodes) is lower than ex- 
pected from random graphs. This arrangement tends to increase the trans- 
mission speed of diseases, energy, matter, and information through networks. 
The largest dista nce in food-web grap hs tends to be small, but the degree of 



clustering varies ()Dunne et al.L 12002 



Pathway proliferation is another large-scale topological characteristic of net- 
works, with implications for energy, matter, and information transmission. It 
is the tendency for the number of pathways in a network to increase geo- 
metrical l y with ou t bound as pa t hway length increases. Patten and colleagues 
()Pattenl . Il985^ : jPatten et al.l Il982 l) first observed this tendency in small, 
well-connected ecosystem model s during the early development of eco s ystern 
network analysis. More recently, Fath ( 1998[ l and Borrett and Patten ( 2003| ) 
showed that the rate of pathway proliferation is variable among networks. 
This is significant because the pathway proliferation rate characterizes how 
quickly the number of indirect pathways increases, and thus, the number of 
pathways available for interactions. Food-web investigations often emphasize 
the shortest pathw ays, assuming that mo s t significaii t inter actions occur over 



these routes (e.g., Caldarelli et al. . 19981 : Post et al. . 2000). However 



previ- 
ous results from Network Environ Analysis, an environmental application and 
extension of economic Input-Output Analysis, indicate that flows over longer 
indirect pathways c an be significant or e ven dominant cons t ituent s of total 
system throughfiow ( Higashi and Patten . 19861 19891 : Patten, 1983| ). which is 
a measure of whole system activity. This result has i mportant impli c ations 
for trophodynamics ( Burns et al. . 1991; Patten et al. . 199^ Whipple . 
and biogeochemica l cycling in ecosvstems (jBorrett et al. . 20061 iFinn . 



1998) 



1980; 



Patten et al. . 1976[ l. Given the possible significance of indirect pathways in 



network models of conservative transport systems like ecosystems, it is critical 
to understand the network characteristics influencing the pathway prolifera- 
tion rate. 



In this paper, we clarify the pathway proliferation concept, describe factors 
influencing the proliferation rate, and characterize pathway proliferation rates 
in 17 large empirical food- webs. In Section 2 we review relevant mathematics 
to build a better understanding of pathway proliferation. In Section 3 we ap- 
ply this understanding to 17 food-web models drawn from the literature. This 
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analysis reveals a type of modularity in some of the food-webs, l ending support 



to the hypothesis that food- webs have a modular structure (lAUesina et al. , 
200,4 iKrause et al.l . l2nn3l : lMavl . ll972l : IPimm and Lawtonill98nl : lYodzij . ll98l ). 



We conclude by summarizing our findings, discussing their relevance for eco- 
logical systems, and suggesting next steps along this research path. 



2 Pathway Proliferation 



Patten et all Il982[ ) intro- 



Although Patten and colleagues (jPattenl . 11985c 
duced pathway proliferation into the ecological literature over two decades ago, 
it is not well understood. In this section, we synthesize mathematical results 
from graph theory and matrix algebra to determine a method for quantifying 
the node-node pathway proliferation rate, and to identify the bounds and ex- 
pected value of the rate. In addition, we uncover the possibility of differing 
rates of pathway proliferation for modules within a network. 

Network models of comp lex systems are mathernatically graphs which can 
be directed or weighted ( Bang- Jensen and Gutin . 2001 : Ponsteinl . Il966[ l. A 
graph G is specified by a set of n nodes and e unoriented edges (0 < e < 
n(n — l)/2 + n), where edges indicate an undirected relationship between two 
nodes. A directed graph (digraph) D is also specified by a set of n nodes, but 
instead of edges it has L oriented arcs or links (0 < L < n^). Digraph structure 
is partially characterized by two connectivity measures, connectance C = Ljr? 
and link density l^ l n-, which are common metrics in the food- web literature 



( Cohen et al.l . 119901 : lMartinea . ll991f ). Edges and links can be assigned weights 



to represent the relationship strength. 



In this paper, we focus on simple unweighted digraphs, where simple implies 
no more than one link from one node to any other. We do this for two reasons. 
First, directed graphs are often appropriate for ecological applications as many 
ecological processes are oriented (e.g., predation and excretion generate energy 
and matter flows from one ecosystem element to another). Second, while the 
mathematics described in this paper may apply to non-simple and weighted 
graphs, our interest here is primarily network structure as it is a necessary 
element of understanding ecosystem organization. 



2.1 Quantifying Pathway Proliferation 



Here, we review relevant definitions and results from graph theory and linear 
algebra. We first show why the dominant eigenvalue of a strongly connected 
digraph is a good measure of the node-node proliferation rate. We conclude 
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by describing how the dominant eigenvalue can distinguish between funda- 
mentally different ecological networks. 



In a directed graph, a pathway is an alternating sequence of nodes and links 
connecting a starting and a terminal node. Pathway length m is the number of 
links in the pathway. For example, in the directed graph D in Figure la there 
is a pathway of length 2 from node 1 to node 3 (e.g., 1 ^ 2 ^ 3). Cycles are 
pathways with the same starting and terminal nodes, and a cycle of length one 
is a self-loop. In our example network, 1— >2— >-3— >lis cycle of length three, 
and 4 — i> 4 is a self-loop. Pathways with self-loops are termed walks, those 
without are paths (Patten, 1985a). Alternatively, we can represent D with its 
associated and isomorphic adjacency matrix A„xn = (O'ij), where a^j = 1 if and 
only if there is a link from j to i (note column to row orientation), otherwise 
ttij = (Figure lb). The number of direct links terminating or starting at a 
node is termed the in- and out-degree, respectively. These are calculated as 



(kr) = E"=i aij and k' 



out 



out^ 

j 



.out 



are 



EILi^ij; where k'" and k° 
1 X n and X 1 vectors respectively. Average in- and out-degrees (< fc*" > or 
< >), and degree distribut i ons Pjk) are ways of c haracterizing network 



struc ture ([Albert and Barabasi . 20021: Newman, 20031 : Ulanowicz and Wolff . 
199lL 



Indirect pathways (m > 1) are enur nerated by rais ing the adjacency ma- 
trix to the m*^ power A™" = (dij)™" (Ponstein, 1966t ). Again, pathway pro- 
liferation is the tendency for the number of pathways in a network to in- 
crease without bound as a function of increasing pathway length (Figure Ic). 
Borrett and Patten (|2003t) approximated whole-system pathway proliferation 
rate as 

^ofl'(EE4r^) 



logiEEaif ^^)) 



where a^-^^ signifies the aij element of A™, but th is agg r egate estimate only 



holds for sufficiently large m. Therefore, following iFathI ( 199a ). here we will 



define pathway proliferation in terms the rate at which a-™''' changes as m 



oo. T his makes pathway proliferation a combinatorial property of A (jSeneta 
19731 ). 



As m increases, the development of A™ is determined by its characteristic 
polynomial, vr(A) = det(AI — A ), where A is the var i able i n the polynomial 
and Inxn is the identity matrix ( Godsil 19931 : Seneta . 1973[ ). The eigenvalues 
(Aj, z = 1, . . . , n) of A are roots of the characteristic polynomial, determined 
as solutions to 7r(A) = 0. The set of eigenvalues IAi > Xp > . . . > Aj > . . . > 
A„} is the spectrum of A, which Cvetkoviac et al. (1980) study to determine 
structural properties of graphs. In undirected graphs, A is symmetric and all 
Aj are real; in directed graphs Aj may be a complex number. 
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Fig. 1. Example digraph, its associated adjacency matrix, and pathway prohfera- 
tion. (a) Digraph D composed of strongly connected components Ki = {1,2,3}, 
K2 = {4}, K3 = {5}, K4 = {6}, K5 = {7}, (b) the adjacency matrix A associated 
with the digraph D (dotted lines denote strongly connected components), and (c) a 
plot demonstrating pathway proliferation from node 3 to 1: the number of pathways 
increases as pathway length m increases. Transient effects created by smaller 
eigenvalues are visible in the first ten pathway lengths, but the rate of pathway 
proliferation has nearly converged to Ai(Ki) = 1.32 by a pathway length of 20. 



In addition to being the roots of the characteristic polynomial, the eigen- 



Wi 



W2 



w. 



IS 



values must satisfy AW = AW, where W^xn = 

a composite matrix of the right eigenvectors Wj, and Anxn = diag(Aj) is 
a matrix with the eigenvalues of A on the principle diagonal and zeros in 
all other positions. Further, the eigenvalues must satisfy VA = VA, where 

is a composite matrix of the left eigenvectors 



V 



Vl 



V2 



V,-. A™ is then determined as 
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^-1 



A'^ = W X A'" X W" 

■■• 




A'" = W 



X 



A™ 



■■■ a: 



X W 



-1 



where W is as before and is the matrix inverse of W, which are the 

transposed complex conjugates of the left eigenvectors V ( Caswelll . 2001 ). If 
we let V* be the i*'' row of W~^, then we can rewrite equation 1 as 



A" 



A^wivt + A^W2V^ + ■ ■ ■ + Arw.v,* + A^Twnv: (2) 

Equation 2 is the pathway generating function of A ( Godsil 1993[ ). and illus- 
trates how the development of A"* depends on the spectrum of A. 



Next, we introduce two classification schemes, one for digraphs and one for 
matrices, because they are required to apply a theorem that will let us de- 
velop a succinct estimate of the pathway proliferation rate. In the first scheme, 
digraphs are classified as one of three types: strongly connected, weakly con- 
nected, and disconnected. A digraph is strongly connected (strong) if it is 
possibl e to reach every node fro r n ever y other over a pathway of unspecified 
length. iBang- Jensen and GutinI ()2nnih define a (sub) digraph with only one 
node as strong, though this is trivial for our purposes. A weakly connected 
(weak) digraph is one in which it is possible to reach any node from any other 
node if we ignore link orientation, but it is impossible when following link 
orientation. Nodes of a weak digraph must have an in- degree or out- degree of 
at least 1. A disconnected graph is one that contains one or more non- adjacent 
strong or weak components. The digraph in Figure la is disconnected. 

The second classification scheme distinguishes three types of adjacency ma- 
trices. A matrix is irreducible if and only if it is associated with a strong di- 
graph, while one associated with a weaklv c onnected or disconnected digraph 
is reducible ( Berman and Plemmonsl . Il979( ) . Irreducible matrices are further 
divided into two classes: primitive and imprimitive. A primitive matrix is an 
irreducible matrix that b ecomes positiv e {aij > 0, for all i,j) when raised to 
a sufficiently large power ( Seneta . 1973[ ). Furthermore, all adjacency matrices 
are non-negative because all aij are greater than or equal to zero. 

Weak and disconnected digraphs are decomposable into a unique set of maximally- 
induced, strong subdigraphs (Kj . i = 1 a, where a < n), th at are termed 

strongly connected components (jBang- Jensen and GutinI . 120011 ). An induced 
subdigraph of D is a subset of nodes in D with all links that both start and 
terminate on the node subset, and a maximally induced subdigraph is the 
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largest one that is strong. This imphes that there is at least one simple cycle 
(no repeated medial nodes) that connects all nodes in a non-trivial strongly 
connected component. Further, adjacency matrices associated with Kj are 
irreducible. For example, the digraph in Figure la contains two connected 
subdigraphs {1,2,3,4} and {5,6,7}. Furthermore, it can be partitioned into 
five strongly connected components Ki = {1,2,3}, K2 = {4}, K3 = {5}, 
K4 = {6}, K5 = {7}, of which only Ki is non-trivial. The adjacency matrix 
associated with each strongly connected component is irreducible; the adja- 
cency matrices associated with Ki and K2 are primitive. 



Given these definitions, the Perron-Frobenius theorem guarantees there is one 
real eigenvalue equal or l arger than all other eigenvalues , A 1 > X j {i = 2, . . . , n) 
in irreducible matrices ( Berman and Plemmonsl . Il979l: ISenetal . Il973l ). In the 
hterature Ai(A) is alternately referred to as the dominant eigenvalue, the 
Perron eigenvalue, and the spectral radius. Next, we illustrate why Ai(A) is a 
good measure of the pathway proliferation rate. 



As shown by Caswelll ( 2001 ). we can divide both sides of equation 2 by Ai to 
obtain. 



A^ 



A5 



a: 



Af 



wivi + -— -W2V2 + --;-W3V3 H h TT;rw„v, 



A? 



A? 



(3) 



If A is primitive and irreducible, then Ai is strictly larger than ||Aj|| for all 
i > 1, where || • || is the norm of • (this is necessary since Aj may be complex). 
Taking the limit of both sides of equation 3 as pathway length increases, we 
find that 

Thus, smaller eigenvalues infiuence pathway proliferation over shorter path 
lengths (Figure Ic), but as path length increases the pathway proliferation 
rate asymptotically be comes Ai. Consequently, A™"*"^/A™' — > Ai as m ^ 00 
f)Hilll . [unpublished ms.i ) . and Ai is the asymptotic rate of pathway proliferation 
in a strongly connected graph with a primitive adjacency matrix. In addition, 
the damping ratio 

WMW 



characterizes the rate of convergence to Ai(A) (|Caswelll . l200lh . 



In strong digraphs with an imprimitive and irreducible adjacency matrix there 
are c < n ei genvalues with the same absolute magnitude, and one or more may 
be complex ( Seneta . 19731 ). The Perron-Frobenius theorem then indicates that 
the common absolute magnitude of the c eigenvalues will be larger than the 
other n — c eigenvalues. In this case, the dominant eigenvalue has a multiplicity 
of c, and as m — » 00 only the c largest eigenvalues will infiuence pathway 
proliferation. Caswelll ( 2001 ) reports that these digraphs generate oscillatory 
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dynamics. 



Given these mathematical results, each strong component Kj of a weak or 
disconnected digraph will have an independent rate of pathway proliferation, 
Ai(Kj) (read "Ai of Kj"). The eigenvalues of a reducible matrix are the union 
set of the eigenvalue s of the adjacency matrices associated with strongly con- 
nected components ( Jain and Krishna . 2003| ). Thus, the maximum dominant 



eigenvalue of the strongly connected components will be the dominant eigen- 
value of the whole digraph. Further, trivial strongly connected components — 
those with only one node — will have a pathway proliferation rate of unity or 
depending on whether or not it has a self-loop. Thus, if a digraph is composed 
of only trivial strongly connected components without self-loops, pathway pro- 
liferation will not occur. This is true of all acyclic digraphs, and suggests that 
we can use the dorainant eigenvalue to detect the presence of cycles in digraphs 
( Jain and Krishna . 20031 ). For this application there are three cases: 



(1) if Ai(A) = 0, then A has no cycles; 

(2) if Ai(A) = 1, then A has at least one cycle and all cycles occur in strongly 
connected components that have only one simple cycle; and 

(3) if Ai(A) > 1, then A has more than one simple cycle. 



Based on an independent development, FathI ()l998| ) interpreted similar results 



as three classes of feedback: 1) no feedback, 2) simple feedback, and 3) cyclic 
feedback in strongly connected networks. Notice that a graph with Ai(A) > 1 
could have a reducible or irreducible adjacency matrix, while the adjacency 
matrix of a graph with Ai(A) = is necessarily reducible with a = n trivial 
strongly connected components. 

Similarly, as m ^ oo we can summarize three possibilities for the dominant 
eigenvalue as a measure of pathway proliferation in digraphs: 

(1) if Ai(A) = 0, then the number of pathways between two nodes declines 
to zero; 

(2) if Ai(A) = 1, then the number of pathways between nodes in a strongly 
connected component remains constant; and 

(3) if Ai(A) > 0, then the number of pathways between nodes in at least 
one strongly connected component (Kj) increases without bound at an 
asymptotic rate equal to Ai(Kj) where max(Ai(Kj)) = Ai(A). 



2.2 Bounds and Expected Values of Pathway Proliferation 



Given that Ai(A) is the asymptotic rate of pathway proliferation in strong 
digraphs, it would be useful to know its theoretical bounds and expected 
value. Again, existing mathematics provides us with some of these answers. 
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Matrix theory bounds the dominant eigenvalue of a non-negative matrix by 
the minimum and maximum column (row) sum, which in the context of di- 
rected graphs is the minimum and maximum out-degree (in-degree), where 
equality holds only if fc*" = k°'^^ (|Berman and Plemmonsl . Il979l : ISenetal . E973.) . 
Thus, max min(A;*"'), min(/c°"*)^ < Ai(A) < min (^max(/c*"'),max(/c°"*)^. In 
a strongly connected digraph with more than one node, all nodes must have 
at least one and a maximum of n incoming and outgoing links. Therefore, 
1 < Ai(A) < n for a strong digraph. As stated previously, a trivial compo- 
nent with no self-loops will have Ai(Kj) = 0, and a complete graph will have 
Ai(A) = n (allowing self-loops). Notice that in the binary matrix A, Ai(A) 
cannot take values between and 1. 

With these bounds, we can now examine the expected value of A(A). In undi- 
rected random graphs G with A = {aij) where aij = aji = 1 with probability 
p {0 < p < 1) and Qij 



with probability (1 — p), Juhasz proved that 



lim 



Ai(A) 



n 



p 



(6) 



( Cvetkoviac and Rowlinson . 1990^. This implies that Ai(A) ~ np in the limit 
of large n. Furthermore, given Ai(A) and n we can determine the approximate 
nur nber of undirected edges in A. A random graph is not necessarily connected, 
but lErdos and Renvil (ll959Lll96nh showed that the fraction of nodes connected 
in a single component increases rapidly when the average link density exceeds 
unity. 

The expected value of Ai (A) is se nsiti ve to the as s umpt ions of random graphs. 
For example, iFarkas et al.l (|2001[ ) and lGoh et al.l (|2001[ ) found that in random 
graphs with a power-law distribution of node degrees rather than the standard 
Poisson distribution, Ai (A) increased with approximately n^^^. Furthermore, 
de Aguiar and Bar- Yam (2QII5) demonstrated that the expected Ai(A) is fur- 
ther modified if the network topology displays a hierarchic modularity. 

Random graphs are well-studied, but properties of random digraphs are less 
well-known. Some characteristics are similar to undirected graphs. For exam- 
ple, the in- and out-degree of random digraphs has a Poisson distribution, 
and when link density is greater than unity the expected size of the l argest 
strongly connected component increases rapidly ( Barbosa et al. . 20031 : Karpl 
19901: lLuczakl . ll99c l). However, we are unaware of results regarding the spectra 
of random directed graphs. Therefore, we numerically verified that Ai(A) ~ np 
remains plausible for random directed graphs by determining the largest eigen- 
value in an ensemble of 99,000 random digraphs (50 from each combination 
of n = {2, 3, . . . , 100} and p = {0.05, 0.10, 0.15, . . . , 1}). Our results indicate 
that Ai(A) ~ nC = L/n, where C = L/n^ is an estimate of p (Figure 2). 
As either L increases or n decreases the residual error decreases. We conclude 
that in random digraphs, as in undirected random graphs, Ai(A) is largely 
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Fig. 2. Dominant eigenvalue in random digraphs, (a) Points show the relationship 
between dominant eigenvalues Ai (A) and link density L/n in 99,000 uniform random 
digraphs where aij = 1 with probability p and Uij = with probability (1 — p) (50 
replicates of each combination of n = {2, 3, . . . , 100} and p = {0.05, 0.10, . . . , 1}). 
The line indicates the expected Ai(A) = L/n relationship, (b) Distribution of resid- 
ual values indicating that as L/n increases it becomes a better predictor of A(A). 

determined by the combination of the number of nodes and number of direct 
connections; pattern of connections has a minor influence. In digraphs with 
a more structured topology — such as those with power-law in-degree or out- 
degree distributions or modularity — we might expect Ai(A) to deviate from 
L/n as it does in undirected graphs, though this remains to be explored. 

In ecological networks where n is the number of species (functional groups, 
etc.) and L is the number of direct transactions, the rate of pathway prolifer- 
ation will be heavily in fluenced by species richn ess and direct liri k abundance. 
However, the results of lFarkas et al.l (|200l[ ) and lGoh et al.l (j200l|) suggest that 
if the degree distributions are skewed, as has been de r nonstrated for some fo od- 
webs f)Dunne et al.1 . 12002^ iMontova and Soi3 . 120021 : IWilliams et aP . 120021) . or 
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the networks contain other types of order such as modularity, then the residual 
error d = |Ai(A) — L/n\ may be larger than expected from random graphs. 

We have presented three key results in this section. First, pathway prolifer- 
ation occurs only if there is more than one cycle or feedback in the graph 
(i.e., it does not occur in acyclic networks). Second, the dominant eigen- 
value of strongly connected components Ai(Kj) is the asymptotic prolifera- 
tion rate for all nodes within that component, and this rate can vary between 
strongly connected components within a network. Further, the largest Ai(Kj) 
of a network is the dominant eigenvalue for the whole network such that 
max(Ai(Kj)) = Ai(A). Third, although topology can be a factor, the prolifer- 
ation rate is heavily influenced by the number of nodes and number of direct 
links. Additionally, networks composed of a non-trivial strongly connected 
components Kj (i = 1, . . . , a) have a form of structural modularity that may 
be functionally significant to the complex systems being modeled. 



3 Modularity and Pathway Proliferation in Food-webs 



In this section we build on our conceptual and mathematical developments 
by applying the theory to determine the pathway proliferation rates in 17 
of the largest empirical food-webs currently available, which are identified 
in Table 1. Five of the food- webs are terrestrial, three are from freshwater 
habitats, and seven represent marine or oceanic habitats. Following accepted 
protocol, original food-webs were modified such that any species or functional 
group with identical predators and prey were grouped as a "trophic spec i es" or 



trophospecies to r e duce naethodo logical bias in the data ( Cohen et al. . IQQOl : 



Pimm et al. . 1991 : Yodzisl Il982 ^. ^998 ). The s e foo d- webs have been the sub- 



ject s everal network analys es (Du nne et al. . l2002a .b. 20041: Williams et al. 



2OO2I : Williams and Martine z. 200_oj) which previously reported their number 
of trophospecies or nodes n, connectance C = L/n'^, the proportion of basal 
species %B {kl"^ = 0), proportion of intermediate species %I {kl'^,k°^^ > 0) 
, proportion of top species %T (/c°"* = 0), and link density L/n. We repro- 
duce this bas i c netw ork information in Table 1 for comparison. In addition, 
(|Dunne et al.l . l2002allbh showed that several have skewed degree distributions 
(i.e., power-law, exponential). 



3.1 Methods 



We first identified and characterized all strongly connected components (Kj), 
including their rates of pathway proliferation and damping ratio, in these 
food-webs. We envisioned three possible outcomes. If food-webs were ade- 
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quately modeled by random digraphs then we would expect each web to have 
one strongly connected component encompassing most if not all of the nodes 
with a single pathway proliferation rate close to link density Ai(A) ~ L/n. 
This seemed unlikely given the known skewed degree distributions and ad- 
ditional evidence that ecological processes c onstr uct non-random topologies 
f Cohen et al. . 199^ Willianas and Martine d . 2000f ). despite arguments to the 
contrary ( Kennv and Loehld . 1991 ). A second pos sibility is based on the ob- 
servation that most early food- webs were acyclic ( Cohen et al. . 199Cl| ). Thus, 
the 17 food-webs in our study, all developed since 1990, might also be acyclic 
digraphs with no non-trivial strongly connected components. In this case, 
pathway proliferation would not occur and Ai(A) = 0. A final possibility is 
that the food-webs would tend to have one or more strongly connected compo- 
nents and multiple pathway proliferation rates. This outcome would support 
the hypothesized modul arity of ecological systerns whi c h is thought to in- 



creas e system stab ility ()Krause et al 
198nl:lYodzij . ll982h . 



20031 : iMavl . Il972l : IPimm and Lawton . 



In Section 2 we hypothesized that the absolute difference between the domi- 
nant eigenvalue and its expected value in random digraphs {L/n) might be a 
useful indicator of the significance of network topology. To assess this hypoth- 
esis, we used Monte Carlo simulations to determine if d = |Ai(A) — L/n\ was 
larger than expected. We had two scales of analysis: whole food-web and non- 
trivial strongly connected components. For both, we constructed 1001 uniform 
random digraphs with n nodes, where each possible link was connected with 
probability p equal to the original network's connectivity {p = C). We as- 
sessed statistical significance by determining the fraction of random digraphs 
in which d was equal or greater than observed in our network of interest, 
Pr{d). Assuming a significance level of a = 0.05, Pr[d) < 0.05 implies d is 
statistically significant. 



When applied to the entire food-web, a significant difference with the null 
model implies that topological factors beyond species and link richness are 
significant in determining the whole system dominant eigenvalue. This could 
be the size or frequency of strongly connected components within the network, 
as suggested by the analysis in Section 2, or perhaps a skewed degree distri- 
bution. If the food- webs had more than one non-trivial component, then we 
expected the deviation to be large. 



When applied to strongly connected components, a significant deviation of d 
also indicates the importance of network topology. However, given that they 
are irreducible by definition, a significant deviation of d within a component 
must indicate the significance of another element of topology, such as the 
degree distribution. 
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3.2 Results 



Food-webs included in this study range in size from 25 to 155 trophospecies 
and 3% to 32% connectance (Table 1). They tend to have a large proportion 
of intermediate species (i.e., those with kj^ > and > 0), although 
the two stream food-webs are notable exceptions. Ten of the 17 food-webs 
examined contained at least one non-trivial strongly connected component; 
six had two. Five of the remaining food-webs had a dominant eigenvalue of 
unity, implying that at least one node contained a self-loop. Our results reveal 
that the majority of these f ood-webs hav e at le ast one directed cycle, contrary 
to earlier food- web theory ( Cohen et al. . 19901 ). 



While the majority of the food-webs have a modular organization that is based 
on strongly connected components, the proportion of species involved in the 
modules is variable. In food-webs that have a non-trivial strongly connected 
component, the proportion of the original nodes involved ranges from 2% in 
the two Ythan Estuary food-webs to 60% in the Caribbean reef model. Notice 
that the definition of a strongly connected component bans nodes that have no 
inputs or no outputs, which excludes basal species {kj^ = 0) and top consumer 
species (A;^"* = 0). Therefore, the total number of species in strongly connected 
components is limited by the number of intermediate species. This may be 
a factor in why the two stream food-webs and the Scotch Broom food-web 
contain no non-trivial components. 

The absolute difference between the dominant eigenvalue of the entire food- 
web and its expected value based on random digraphs of the appropriate 
size and connectance {d = |Ai(A) — L/n\) ranged from a minimum of 0.5 
for the El Verde rainforest to a maximum of 12.9 for the NE US shelf food- 
web (Table 1). In all cases, this difference was significantly different from the 
random digraph null model, indicating that topology is a significant factor 
in determining Ai(A) . This result is consistent with the presence of one or 
more small non-trivial components and acyclic digraphs. Our results provide 
another line of evidence suggesting that the ecological processes that create 
food-webs lead to more ordered network topologies; random digraphs are not 
good models for these systems. 

Inspection of the strongly connected components reveals a diversity of topolo- 
gies as shown in Table 2. The largest strongly connected component, with 
seventy trophospecies occurs in the El Verde rainforest model; although it is 
the least well connected (13%), it still has the largest rate of pathway pro- 
liferation (Ai(A) = 10.25). In contrast, ten of the sixteen strongly connected 
components have four or fewer species. Five of the strongly connected compo- 
nents only contain two trophospecies, requiring a single simple cycle of path 
length 2 (e.g., j — > i — > j). While the two strongly connected components in 
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the Coachella Valley, Skipworth pond, and Benguela are about the same size, 
one of the two components in Little Rock Lake, Caribbean reef, and NE US 
shelf is substantially larger than the other. Table 3 lists the trophospecies in 
the two strongly connected components of the Coachella Valley. 

The dominant eigenvalues of all strongly connected components have a mul- 
tiplicity of one, so the adjacency matrices associated with the component 
subdigraphs are primitive. Therefore, the dominant eigenvalues represent the 
strongly connected component asymptotic rates of pathway proliferation. These 
range from 1.62 in strongly connected components of the two Ythan Estu- 
ary food-webs to 10.25 in the large El Verde rainforest strongly connected 
component and generally increase with link density as would be expected in 
random digraphs. However, half of the strongly connected components have 
a statistically significant difference between the dominant eigenvalue and link 
density including Coachella Valley (Ki and K2), El Verde Rainforest (Ki), 
Little Rock Lake (Ki), Benguela (K2), Caribbean Reef (K2), and NE US shelf 
(Ki and K2). The topological arrangement of species and links in these three 
strongly connected components influences their rate of pathway proliferation; 
the others are largely determined by their species and link richness. 



The damping ratio defined in equation 5 is an index of the speed of conver- 
gence to the asymptotic rate of pathway proliferation; A larger ratio indicates 
faster convergence. Five of the strongly connected components are completely 
connected. They have a pathway proliferation rate equal to their trophospecies 
richness and an undefined damping ratio because their second eigenvalues are 
zero. In these cases the asymptotic rate of pathway proliferation is achieved 
instantaneously. The other damping ratios range from 1.69 in K2 of the NE US 
Shelf to 4.41 in Ki of the same food-web. Transient dynamics of the pathway 
proliferation rate, determined by the smaller eigenvalues, are more influen- 
tial in NE US shelf (K2). Its pathway proliferation rate does not converge 
until a pathway length of nearly twenty-two, while in NE US shelf (Ki) the 
proliferation rate converges by a pathway length of eight. 



In summary, the majority of the food-webs we examined contained at least one 
non-trivial strongly connected component. Six food- webs had two non-trivial 
strongly connected components; none had more than two. The proportion of 
species involved in a strongly connected component ranged from 2% to 60%. 
In all cases, the difference between the dominant eigenvalue of the food-web 
and the expected value (L/n) in a random network was significant. This differ- 
ence occurs because the topology of food- webs is non-random; thermodynamic 
processes and species ch aracteristics (e.g., the spe cies niche) combine to form 
non-random structures ( Chase and Leibold . 2003| ). Within the strongly con- 
nected components, the rate of pathway proliferation ranged from 1.62 to 10.25 
and half were indistinguishable from random graphs based on the expected 
rate of pathway proliferation. 
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4 Discussion 



As with any analysis of network models that reveals previously hidden struc- 
tural patterns, we are left with two questions. First, what, if any, significance 
do these patterns hold for our systems of interest? Are strongly connected 
components and pathway proliferation simply another pretty pattern, another 
network or food- web statistic to report, or do they impart some functional sig- 
nificance? The second question cannot be divorced from the first; what system 
processes might create these structural patterns? Are there ecological processes 
or forces that might lead to the development of these structures? These are 
not easy questions to answer, but in this section we attempt to address them 
for the presence of strongly connected components and pathway proliferation 
in ecological networks. 



Strongly connected components introduce a form of modularity into network 
models, where modularity is defined as a hierarchical system subdivision into 
more or less interacting subsystems. Several types of modularity have been 
proposed in ecological systems. Building on earlier ideas in general systems 



theory that linked hierarchical organization to system stability (e.g. ISimon 



19621 ). Mav ( 1973| ) hypothesized that ecosystems have modular structures. He 



found that the Lyapunov stability of randomly assembled ecosystems tended to 
be greater when the species were partitione d into blocks of interacting species 
with few if any connections to other blocks. Pimm ( 1979[ l termed these blocks 
of species "compartments" , stating that they are "...characterized by strong in- 
teractioi is within compartments, b ut weak interactions among compartments" 
(p. 145). IPimm and Lawtonl (|l98(i concluded from a study of binary empirical 
food- webs that, while there was evidence species were grouped into subsystems 
largely by habitat, compartmentalization as defined by Pimm was an uncom- 
mon phenomenon. They noted, however, that a complete test of the hypothesis 
would require knowledge of the strength of interactions, which was absent in 
their food-web models. It is also possible that they were unable to identify 
compartmentalization in their food -webs because their models exclude detri- 
tal recycling. However, iNeutel et al, (.2002.) showed that even without detrital 
recycling, long, weak links enhanced system stability. 



Yodzis (1982) remarked that modula r organization was an old idea in ecology, 
citing the guild concept (|Rootlll967l ) as an example. He applied the dominant 
clique idea from graph theory to identify another type of modularity in food- 
webs. He defined "a clique as a set of species in a given ecosystem with the 
property that every pair in the set has some food resource in common, and 
... a dominant cliqu e as a clique which is contained in no other clique" (p. 
552). More recently, iKrause et aP (|2003h used a methodology developed to 
identify cohesive subgroups in social network analysis to classify another type 
of modularity in food-webs. They demonstrated that this type of organization 
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increased system sta bility to specie s delet ion by localizing the effect within a 
module. In addition, Allesina et al. ( 20051 ) found four or more modules based 
on strongly connected components in 17 network models of carbon flux. The 
four modules usually grouped into four types: inputs, outputs, dissipation, and 
species and nutrient pools. Their work is the most similar to that presented 
in this paper, but there are two important differences. First, we examined 
different types of ecological networks. They examined ecosystem flow networks, 
while we restricted our analysis to food-webs. Second, we considered inputs, 
outputs, and dissipation to be external to the system and therefore we could 
not identify these as separate modules. 



The dominant ecological hypothesis is that food-web or ecosystem modularity 
increases overall system stability by localizing interactions within modules. 
Given the static, binary, presence-absence information of food-webs in our 
study we were unable to meaningfully test this hypothesis; stability is inher- 
ently a dynamic concept. Known issues with food- web model construction fur- 
ther m ake this hypothesis difficult to resolve ( Cohen et al. . 19931: Polls and Strong . 
199fit ). Empirical food-web models usually do not indicate interaction strength 
or the temporal and spatial variation of the interactions, as these details are 
expensive to acquire. 



Despite the challenge of assessing their effect on system stability, strongly 
connected components in ecological networks appear to be important func- 
tional elements of system organization and provide new insights about species 
participating in them. By definition (Section 2), there is minimally one sim- 
ple cycle that encircles all nodes in the strongly connected component. This 
provides at least one channel for cybernetic feedback (positive or negative) 



to propagate among species in the module ( DeAngelis et al. . 19861: Patten 



I959I: Patten and Odum . 1981 ). Furthermore, in food- webs it is reasonable 
to assume that each predatory species directly benefits by its consumption 
of prey. This establishes an indirect mutualism that spans the strongly con- 
nected component, and provides the necessary conditions for the strongly con- 
nected components of food-webs to fu nction as autocatalytic cycles — systems 



that catalyze th eir own production ([Mavnard Smith and Szathmarvl . 119951 : 
UlanowicZ 1997). Autocatalytic cycles are an essential element of metabolism 
in c hemical and living systems and r nay h a ve played a role in the origin of 
l ife (iMavnard Smith and Szathmarv . W9^. Mavnard Smith and Szathmarv 
describe autocatalytic cycles as a force for cooperation among the mem- 
ber species and efficient information integrators. 



Ulanowic3 (1997) identifies several emergent properties autocatalytic cycles 
may possess, including centripetality, persistence, and autonomy. Centripetal- 
ity is the tendency of the cycle to attract more of the energy-matter flux of the 
system. If any member species becomes more efficient at using its resources or 
better able to acquire new resources such that its population increases, this 
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positive change tends to cascade through the module, collectively benefiting 
the populations of all species involved. Autocatalytic cycles tend to persist 
because their general form can be maintained in a system despite fiuctuations 
in the interaction strengths and possible element replacement. In food-webs 
this implies that when a trophically similar species is introduced to the sys- 
tem, if it is more efficient or in some way ecologically more competitive, it 
may wholly replace an existing species in the autocatalytic cycle, but the cy- 
cle remains. Finally, autocatalytic cycles can establish a degree of autonomy 
because species in the cycle can actively infiuence at least a portion of their 
input environment. In this sense, species in strongly connected compo nents 



of food- webs are in volved in ecosystern engineering d Jones et al. . 199?I)_and 
niche constr uction ( Laland et al. . 19991 : Odling-Smee et al. 1 120031 ). luilnowiczl 



f 19861 1997 ) further argues that the autocatalytic nature of cycles in ecosys 
tems makes them a principal agent in ecosystem growth and development. 



Mavnard Smith and Szathmarv (|199,^ remarked that autocatalytic cycles are 
sensitive to cheaters or parasites that feed off the strongly connected compo- 
nent without participating in the cycle. Top predators feeding on species in 
a strongly connected component or a downstream strongly connected compo- 
nent might function as parasites in this sense. Perhaps this is why strongly 
connected components do not occur in all the food-webs analyzed. This and 
the tendency for centripetality may explain why there are fewer than three 
strongly connected components in these food webs. However, we are unable to 
assess these possibilities with these data because the differences may simply 
reflect disparities in food-web modeling decisions. 



Pathway proliferation rates of strongly connected components provide addi- 
tional information about these modules. Each additional link in a strongly 
connected component beyond those that form the deflning cycle introduces 
another embedded simple cycle. This lowers the maximal distance between 
nodes in the module, increases the potential pathways for energy-matter flux, 
tends to increase the rate of pathway proliferation, and leads to the unbounded 
growth of pathways as length increases. In some cases, the rate of pathway 
proliferation will not increase as expected with the number of links. For exam- 
ple, half of the strongly connected components identifled in our food-webs had 
pathway proliferation rates that were signiflcantly different from the expected 
rate based on random graphs with a Poisson degree distribution. This sug- 
gests that module topology differs from what we would expect from a random 
generating process. 

As mentioned earlier, we expected the ecological processes forming food-webs 
to generate non-random structures. Species characteristics such as metabolic 
requirements, food preferences, capture ability and handling time, as well as 
other niche requirements and natural history constraints should directly influ- 
ence the choices of "who eats whom" and how much. In addition, the emergent 
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properties of autocatalytic cycles and ecosystems more generally may provide 
whole-system constraints. 

Perhaps the interesting question is not, why the eight strongly connected com- 
ponents did not match the random expectation, but why the other half did? 
Notice that the strongly connected components with non-random topologies 
were the largest modules, while the eight strongly connected components with 
topologies indistinguishable from random digraphs involved only two or three 
trophospecies. The universe of possible topologies is much smaller in these 
small and well-connected (0.75 < C < 1) modules, making the ecologically 
created topologies reflected in the food- webs more likely. Five of these modules 
were completely connected, generating only one possible arrangement. The 
eight strongly connected components with apparently non-random topologies 
were less well connected (0.13 < C < 0.78), generating a much larger uni- 
verse of possible topologies. In six cases the pathway proliferation rate was 
significantly less than expected, but in the large strongly connected compo- 
nent of the El Verde rainforest and the smaller module of the Caribbean Reef, 
the pathway proliferation rate was more than expected. At this point, we are 
unable to provide a satisfactory explanation for these differences. 

In our discussion thus far, we have been interchanging food-webs and ecosys- 
tems. It is important to recognize, however, that food-webs are a subset of 
a broader class of ecosystem models of energy-matter flow. Food-wcbs are 
generally defined by the relation "who eats whom" that is one process gener- 
ating energy-matter flux, while ecosystem flow-storage models typically trace 
a conserved flow unit (e.g., energy, nitrogen, phosphorus) through the sys- 
tem, regardless of the process producing the flow. Thus, non-trophic ecologi- 
cal processes such as excretion and death are captured in flow-storage models, 
revealing a different picture of ecosystem organization with additional cycling. 

Pathway proliferation influences the development and significance of indirect 
flows in ecosystem flow-storage models. Indirect flows are an important as- 
pect of the ecological signiflcance of the strongly connected components, so 
here we take a closer look. Indirect flows are derived from two fundamentally 
distinct types of pathways: chains (e.g., 5 — > 6 — > 7 in Figure la) and cycles 
(e.g., 1 — > 3 1 in Figure la). Indirect flows in chains are limited by transfer 
efficiencies and chain length. In cycles, the number and length of pathways 
arc unlimited such that indirect flows are only limited by transfer efficiencies 
reflecting energy-matter dissipation and export. As ecosystems are open ther- 
modynamic systems, shorter indirect pathways individually will tend to carry 
larger indirect flows than longer indirect pathways. A faster rate of pathway 
proliferation Ai(A) implies that there will be more shorter indirect pathways, 
increasing the possibility that the magnitude of indirect flows will surpass that 
of direct flows. Thus, within a strongly connected component the large num- 
ber of indirect pathways will tend to carry a large fraction of the flow between 
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species (nodes). 



More generally, Ai(A) indicates the potential for direct and indirect energy, 
matter, and information transmission between compartments in a strongly 
connected component. Realized transmission rates are dependent on the re- 
alized use of each pathway. Previous ecosystem network a nalyses reveal some of 

the system- l evel consequences of differential pathway us e (iFathLbnnilFath and Pattenl . 
" 1999I : iHigashi and Pattenl . Il989l : IPattenl . Il985bl : lUlanowiczL Il98fih . but 
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there is much left to learn about this subject. The interplay of this potential 
and realized network structure is an interesting, important topic for under- 
standing the organization and transformation of complex adaptive systems 
like ecosystems. 



We conclude that the strongly connected components and pathway prolifera- 
tion are ecologically relevant phenomena because they provide novel insights 
about the system of interest. Without knowing the strength of interactions 
or energy-matter flux rates, the presence of these structural features suggests 
groups of species functionally integrated by indirect effects mediated by auto- 
catalytic cycles. They portend the possibility of integral species relationships 
that are shifted toward more positive associations and the possibility of the 
dominance of indirect flows. In some cases, apparent negative interactions 
such as predation or competition may become more positive through indirect 
interactions mediated by the autocatalytic cycles of the strongly connected 
components. 
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Table 1. Topological properties of 17 empirical food-webs. 
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Taxa refers to the original number of species; n is the number of nodes or trophospecies; C = L/n"^ is connectance; %B, %I 
and %T are the proportions of basal (indegree = 0), intermediate (indegree and outdegree > 0), and top (outdegree = 0) 
trophospecies; L/n is link density; Ai(A) is the dominant eigenvalue of the entire foodweb; d = |Ai(A) — L/n\] Pr{d) is the 
fraction of an ensemble of 1001 random digraphs in which d is greater than or equal to that observed in K,; * indicates 
statistically significant Pr{d) at a = 0.05; #K is the number of non-trivial strongly connected components; %K is the percent 
of species in a non-trivial strongly connected compo r ient. t marks topological properties previously reported for these food-webs 
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Table 2. Topological properties of strongly connected components in 17 empirical food-webs 
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Kj indicates the non-trivial strongly connected component number; n is the number of nodes (trophospecies), L is the number of links; 
C = Ljn^ is connectance, Ljn is the link density, Ai(Kj) is the dominant eigenvalue of Kj, mult(Ai) is the multiplicity of the dominant 
eigenvalue, p = Ai(Kj)/||A2(Kj)|| is the damping ratio (- indicates p is undefined because |[A2(Kj)|| = 0), d = |Ai(A) — L/n\ is the 
absolute difference between the dominant eigenvalue and link density, Pr{d) is the fraction of an ensemble of 1001 random digraphs in 
which d is greater than or equal to that observed in Kj, and * indicates statistically significant Pr{d) at a = 0.05. 



Table 3 

Trophospecies in the two non-trivial strongly connected components of the Coachella 
Valley food-web. 



Ki 

primarily herbivorous mammals and birds 
small omnivorous mammals and birds 
primarily carnivorous lizards 
primarily carnivorous snakes 
large primarily predacious birds 
large primarily predacious mammals 



K2 

small arthropod predators 
medium arthropod predators 
large arthropod predators 
facultative arthropod predators 
life-history arthropod omnivore 
spider parasitoids 
primary parasitoids 
hyperparisitoids 
predacious mammals and birds 
arthropodivorous snakes 
primarily arthropodivorous lizards 
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